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Abstract 

Nonlinear coupling between 3-minute oscillations and Alfven waves in the solar 
lower atmosphere is studied. 3-minute oscillations are considered as acoustic waves 
trapped in a chromospheric cavity and oscillating along transversally inhomogeneous 
vertical magnetic field. It is shown that under the action of the oscillations the tem- 
poral dynamics of Alfven waves is governed by Mathieu equation. Consequently, the 
harmonics of Alfven waves with twice period and wavelength of 3-minute oscillations 
grow exponentially in time near the layer where the sound and Alfven speeds equal 
i.e. c s ~ va- Thus the 3-minute oscillations are resonantly absorbed by pure Alfven 
waves near this resonant layer. The resonant Alfven waves may penetrate into the 
solar corona taking energy from the chromosphere. Therefore the layer c s ~ va may 
play a role of energy channel for otherwise trapped acoustic oscillations. 
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1 Introduction 



Waves play an important role in the dynamics of the solar atmosphere. They 
are believed as carriers of photospheric energy into the corona leading to 
plasma heating. Waves are intensively observed in the solar chromosphere and 
corona by SOHO (Solar and Heliospheric Observatory) and TRACE (Transi- 
tion Region and Coronal Explorer) (Nakariakov & Verwichte 2005). Dominant 
intensity oscillations in a chromospheric network have a period of ~3-min. The 

Email addresses: dato.k@genao.org, temury@genao.org (D. Kuridze and T.V. 
Zaqarashvili). 



Preprint submitted to Elsevier 



5 February 2008 



3-min oscillations are believed to be either standing waves in a chromospheric 
cavity (Leibacher & Stein 1981) or propagating waves, which arise as the re- 
sponse of the atmosphere to a general disturbance set at its base (Fleck & 
Schmitz 1991). Coupling of these oscillations to the Alfven waves is a cru- 
cial process in both the Alfven waves may easily pass through the 
transition region taking the energy into the corona. 

On the other hand, recent numerical (Rosenthal et al., 2002), analytical (Za- 
qarashvili & Roberts, 2006) and observational (Muglach et al., 2005) studies 
point out the importance of (3 ~ 1 region of the solar atmosphere in the wave 
coupling phenomena. 

Here we consider the 3-min oscillations as the standing acoustic waves oscillat- 
ing along an uniform vertical magnetic field and study their weakly nonlinear 
coupling to Alfven waves near the chromospheric (3 ~ 1 regions. 



2 The wave coupling 

We use Cartesian coordinate system (x,y,z), where spatially inhomogeneous 
(along the x axis) magnetic field is directed along the z axis, i.e. Bq = 
(0, 0, B z (x)). Unperturbed plasma pressure po(x) and density po(x) are also 
assumed to vary along x. In the equilibrium, magnetic and hydrodynamic pres- 
sures satisfy the transverse balance, i.e. Po(x) + B^(x)/8ir = const. Plasma 



(3 is defined as (3 = 8ttp (x)/B^(x) = 2c\j^v 2 A [x\ where c s = J^po/po and 



va(x) = B z /y/4:7ipQ are the sound and Alfven speeds respectively and 7 is the 
ratio of specific heats. For simplicity, temperature is suggested to be homoge- 
neous so that the sound speed does not depend on the x coordinate. 

We consider wave propagation along the z axis (thus along the magnetic field) 
and wave polarisation in yz plane. Then in this case only sound and linearly po- 
larised Alfven waves arise. The velocity component of sound wave is polarised 
along the z axis and the velocity component of Alfven wave is polarised along 
the y axis. Then the ideal MHD equations can be written as 






(2) 



where p = Po + pi and p = po + Pi denote the total (unperturbed plus per- 
turbed) hydrodynamic pressure and density, u y and u z are velocity perturba- 
tions (of the Alfven and sound waves, respectively), b y is the perturbation in 
the magnetic field. Note, that the x coordinate stands as a parameter in these 
equations. 

Eqs. (l)-(5) describe the fully nonlinear behavior of sound and linearly polar- 
ized Alfven waves propagating along the magnetic field. However the sound 
waves may be trapped in a chromospheric cavity between the photosphere and 
transition region leading to standing patterns. Therefore here we consider the 
sound waves oscillating along the z axis and bounded at z = and z = I 
points. Thus there are the standing patterns 

u z = v(t) sin(k n z), p 1 = p(t) cos(k n z), (6) 

where k n is the wavenumber of sound wave such that k n l = 2irl/X n = nir, so 

— n /2, where n — 1, 2... denotes the order of corresponding harmonics. 

It is natural that almost whole oscillation energy in bounded systems is stored 
in a fundamental harmonic. Therefore here we consider the first (n = 1) har- 
monic of acoustic oscillations, however the same can be applied to harmonics 
with arbitrary n. Recently, Zaqarashvili and Roberts (2006) found that the 
harmonics of acoustic and Alfven waves are coupled when the wave length 
of acoustic wave is half of Alfven wave one. Therefore let express the Alfven 
wave components as 

u y = u(t) sin(kAz), b y = b(t) cos(A;^), (7) 

where k^ is the wavenumber of the Alfven waves and the condition k± = 2fc^ 
is satisfied. 

Then the substitution of expressions (6)-(7) into Eqs. (l)-(5) and averaging 
with z over the distance (0, 1) leads to 

db 

— = k A B u + k A vb, (8) 



du 
dt 



kAB 
47rp 



b — k A uv, 



(9) 



dp 

m = ~ pohv - 



'ID 



Here Eqs. (10)-(11) describe the time evolution of acoustic oscillation forced by 
ponderomotive force of Alfven waves, while Eqs. (8)-(9) governs the dynamics 
of Alfven waves forced by acoustic waves in a parametric way. It must be 
noted that the coupling between sound and Alfven waves at (3 ~ 1 has been 
recently studied by Zaqarashvili & Roberts (2006). They consider the general 
case of propagating waves and show an alternate energy exchange between the 
waves during the propagation. Here we consider the coupling between standing 
patterns of the waves, which is the particular case of their results. 



Substitution of u from Eq. (8) into Eq. (9) and neglecting of all third order 
terms leads to the second order differential equation for Alfven waves 



d 2 b 



Po 



b = 0. 



(12) 



This equation reflects the parametric influence of standing acoustic wave due 
to the density variation. Then the particular time dependence of density per- 
turbation determines the type of equation and consequently its solutions. If we 
consider the initial amplitude of Alfven waves smaller than the amplitude of 
acoustic waves, then the term with b 2 in Eq. (10) can be neglected. Physically 
it means that the back reaction of Alfven waves due to the ponderomotive 
force is small. Then the solution of Eqs. (10)-(11) is just harmonic function 
of time p = apocos(uit), where uj\ is the frequency of the first harmonic of 
standing acoustic wave and a > is the relative amplitude. Here we consider 
the small amplitude acoustic waves a < 1, so the nonlinear steepening due to 
the generation of higher harmonics is negligible. Then the substitution of this 
expression into Eq. (12) leads to Mathieu equation 

2 b 

— + k\v\ [1 - 2a cos(wit)] b = 0. (13) 



The solution of this equation with the frequency of o>i/2 has an exponentially 
growing character, thus the main resonant solution occurs when 

u A = v A k A = y , (14) 



where u A is the frequency of Alfven waves. Since k A = ki/2, the resonance 

takes place when v A = c s . Since the Alfven speed v A (x) = B (x)/ ^47ip (x) is 
a function of the x coordinate, then this relation will be satisfied at particular 
locations along the x axis. Therefore the acoustic oscillations will be resonantly 
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Fig. 1. Numerical simulations of wave conversion in c s = va region. Upper panels 
show the density and velocity components of standing acoustic oscillations, while 
lower panels show the magnetic field and velocity component of Alfven waves. Rel- 
ative amplitude of acoustic oscillations is 0.1. Alfven waves with twice the period 
of acoustic oscillations show rapid exponential increase in time. 

transformed into Alfven waves near this layer. We call this region swing layer 
(see similar consideration in Shergelashvili et al. 2005). 

Under the resonant condition (14) the solution of Eq. (13) is 

7 / \ 7 \ r u)\ . u\ i 

b{t) = boexp\ — - — tj cos— t — sin— r , (15) 



where b = 6(0). 

The solution (15) has a resonant character within the frequency interval 
\ua — ui/2\ < \aui/2\. This expression can be rewritten as \va/c s — l\ < \a\. 
Thus the thickness of the resonant layer depends on the acoustic wave am- 
plitude. Therefore the acoustic oscillations are converted into Alfven waves 
not only at the surface va = c s but in the region where c s (1 — a) < Va < 
c s (1 + a). Thus the resonant layer can be significantly wider for stronger am- 
plitude acoustic oscillations. Note, that the resonant Alfven waves expressed 
by Eqs. (7) and (15) are standing patterns with the velocity node at the bot- 
tom boundary (z = 0) and antinode at the top boundary {z — I); the wave 
length of Alfven waves is twice than that of the acoustic oscillations due to 
the condition fc^ = ki/2. Therefore the oscillation of magnetic field lines at 
the upper boundary may excite the waves in an external plasma, which carry 
energy away. 
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Fig. 2. The same as on Fig.l, but for c s = region. Amplitude of Alfven waves 

(lower panels) shows no increase in time, so there is no wave coupling in this region. 

It must be noted that the Alfven waves may undergo a phase mixing due to 
the x dependence of the Alfven speed (Tsiklauri & Nakariakov 2002). However 
the aim of this paper is to show the energy conversion from 3-min oscillations 
into Alfven waves only, therefore we do not consider the effect of phase mixing 
here. 

Numerical solutions of Eqs. (8)- (11) (here the back reaction of Alfven waves, 
i. e. the second term of right hand side in Eq. (10), is again neglected) are pre- 
sented on Fig. 1-2. Figure 1 shows the wave dynamics in c s = va region and we 
see the rapid growth of Alfven waves with twice the period of acoustic oscilla- 
tions. However Fig. 2 shows the wave dynamics away from the resonant layer 
(in c s = 0.7va region) and we see no energy exchange between the waves. Thus 
there is a good agreement between numerical results and analytical solutions. 

It must be mentioned that the equilibrium used in this paper is simplified 
as gravitational stratification, which is important in the solar chromosphere, 
is ignored. The stratification leads to Klein-Gordon equation for propagating 
waves with a cut-off for wave frequencies (Roberts, 2004). The 3-min oscil- 
lations and Alfven waves are above the cut-off, so they may propagate in 
the chromosphere. Unfortunately, the stratification greatly complicates math- 
ematical description of the non-linear coupling. On the other hand, plasma (3 
can be constant along a vertical magnetic tube even in the case of stratifica- 
tion (Roberts, 2004); this is the case when sound and Alfven speeds vary with 
height but their ratio remains constant. Therefore if the waves are coupled at 
/3 ~ 1 region in an unstratified atmosphere, then the same can be expected in 
a stratified medium. However the wave coupling in the stratified atmosphere 
needs further detailed study, but is not the scope of present paper. Our goal is 
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just to show that the 3-min oscillations may transfer energy into incompress- 
ible Alfven waves in f3 ~ 1 region, which was recently observed by Muglach et 
al. (2005). 



3 Conclusions 

Here we show that acoustic oscillations trapped in transversally inhomoge- 
neous medium can be resonantly absorbed by Alfven waves near the layer 
of Va ~ c s . The spatial width of the layer depends on amplitude of acoustic 
oscillations and can be significantly wider for strong amplitude oscillations. 

We consider the observed 3-min oscillations as the standing fundamental har- 
monic of acoustic waves trapped in the solar chromospheric cavity with vertical 
magnetic field and show that their nonlinear coupling to Alfven waves may 
take place near (3 ~ 1 layer. The coupling may explain the recent observational 
evidence of compressible wave absorption near f3 ~ 1 region of solar lower at- 
mosphere (Muglach et al., 2005). The amplified Alfven waves with the period 
of ~ 6 min may carry energy into the corona. There they may deposit the 
energy back to density perturbations leading to observed intensity variations 
in coronal spectral lines. 

Thus the layer of f3 ~ 1 may play a role of energy channelioi otherwise trapped 
acoustic oscillations guiding the photospheric energy into the solar corona. 
Therefore the process of wave coupling can be of importance for coronal heat- 
ing problem, but requires further study especially for a stratified atmosphere. 
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